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Abstract 

Research  was  performed  to  verify  previously  obtained 
friction  coefficients  for  compressible,  laminar  and  turbulent 
flow  In  a  simulated  heat  pipe  for  Intermediate  mass  Injection 
and  suction  rates.  The  research  was  divided  Into  two  main 
parts:  an  experimental  study  and  a  numerical  study. 

In  the  experimental  study,  a  simulated  heat  pipe  con¬ 
sisting  of  a  long  porous  tube  with  air  Injected  Into  one  end 
to  simulate  the  evaporator  and  extracted  from  the  other  end 
to  simulate  the  condenser  was  Investigated.  Axial  static 
pressure  variations  for  various  air  supply  mass  flow  rates 
were  measured.  It  ms  found  that  the  friction  In  the  pipe 
decreased  the  amount  of  pressure  recovery  at  the  end  of  the 
pipe  and  that  this  decrease  became  even  greater  with  Increased 
supply  mass  flow  rate. 

The  pressure  distributions  obtained  In  the  experimental 
study  were  used  as  Input  for  a  numerical  simulation  In  the 
numerical  portion  of  this  research.  A  steady,  one-dimensional 
computer  code  was  developed  to  simulate  the  porous  pipe  sys¬ 
tem  used  In  the  experimental  study  in  order  to  calculate  Mach 
numbers,  axial  Reynolds  numbers,  and  friction  coefficients 
along  the  pipe»v  Newton's  second  law,  continuity,  and  Shapiro's 
method  of  influence  coefficients  were  used  in  the  computer 
simulation  to  calculate  the  flow  properties. '^The  results 


obtained  for  the  friction  coefficients  in  the  porous  pipe 
system  studied  did  not  correlate  with  the  values  predicted 
by  Bowman  *  s  expressions.  They  dl4r  however,  exhlbit<,the 
general  trend  expected  for  the  friction  coefficient  -  axial 
Reynolds  number  product  with  Mach  number . 


COMPRESSIBLE  FRICTION  COEFFICIENTS  IN  A 


SIMULATED  HEAT  PIPE 


Lx.  Introduction 

The  purpose  of  this  introductory  chapter  is  to  acquaint 
the  reader  with  the  history  and  operation  of  the  heat  pipe 
and  the  goals  of  this  research.  The  chapter  is  divided  into 
three  sections.  The  first  section  provides  a  brief  discussion 
of  the  development  and  operation  of  the  heat  pipe,  while  in  the 
second  section  the  pertinent  literature  is  reviewed.  In  the 
last  section,  the  objectives  of  this  study  are  outlined. 

EeyclcLprocQt  and.  Operation  a£  the  Heat  Pipe 

The  concept  of  the  heat  pipe  was  first  conceived  by 
R.S.  Gaugler  in  1942  (4:1).  The  first  serious  development  work 
took  place  many  years  later  by  G.M.  Grover  and  his  associates 
at  the  Los  Alamos  Scientific  Research  Laboratory  in  1964  (4:1). 
Grover  showed  that  the  heat  pipe  is  an  effective  high-perform¬ 
ance  heat  transmitter,  developed  some  applications  for  it,  and 
gave  it  its  present  name.  Following  Grover's  work.  Interest 
in  the  heat  pipe  quickly  grew  for  both  terrestrial  and  space 
applications,  resulting  in  numerous  papers  and  several  books 
being  written  on  the  subject.  Today  heat  pipes  have  numerous 
applications.  Some  of  these  applications  Include  cooling 
nuclear  reactors,  infrared  sensors,  laser  systems,  generators. 
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electronic  and  electrical  equipment,  and  turbine  blades.  Heat 
pipes  have  also  been  used  for  cyrogenlc  eye  surgery.  In 
spacecraft  for  temperature  control  of  spacesults  and  vehicles, 
and  for  heat  collection  from  solar  energy,  geothermal  energy, 
and  exhaust  gases. 

The  operation  of  the  heat  pipe  Is  similar  to  that  of 
the  thermal  syphon,  shown  in  Figure  1.1.  The  thermal  syphon 
consists  of  a  sealed  tube,  from  which  the  air  has  been 
evacuated,  with  a  small  quantity  of  liquid  placed  inside. 

When  the  end  of  the  tube  is  heated,  it  causes  the  water  to 
vaporize  and  move  to  the  cooler  end  of  the  tube  where  It 
condenses  on  the  walls.  The  condensate  then  returns  to  the 
other  end  of  the  tube  by  means  of  gravity.  The  thermal  syphon 
is  limited  in  its  orientation  since  the  evaporator  region  must 
be  at  the  lowest  end.  The  heat  pipe,  however,  consists  of  a 
closed  tube  with  a  porous  capillary  wick  fixed  to  the  inside 
surface,  in  its  conventional  form  (Figure  1.2).  The  wick 
contains  the  liquid  phase  of  the  working  fluid  while  the 
remaining  volume  of  the  tube,  known  as  the  vapor  region, 
contains  the  vapor  phase  of  the  working  fluid.  When  heat  is 
applied  to  evaporator  end,  the  fluid  in  that  section  of  the 
pipe  is  vaporized.  This  creates  a  pressure  difference  down 
the  pipe  and  causes  the  vapor  to  flow  from  the  evaporator  end 
to  the  condenser  end  where  the  vapor  condenses  and  releases 
its  latent  heat  of  vaporization  to  a  heat  sink.  The  evapora¬ 
tion  process  depletes  the  liquid  in  the  evaporator  section 
and  causes  the  liquid-vapor  interface  to  move  into  the  wick 
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have  been  written  on  heat  pipe  theory  and  operation.  This 
study  Is  concerned  primarily  with  the  effect  of  mass  injection 
and  extraction  on  the  friction  coefficient  at  the  tube  wall. 


Several  articles  have  been  written  on  this  subject,  but  only 
a  few  of  these  will  be  discussed  in  this  section. 

Kinney  (5)  studied  the  effect  of  wall  mass  transfer  on 
the  frictional  and  heat  transfer  characteristics  of  fully 
developed  laminar  flow  in  circular  porous  tubes.  His  results 
are  based  on  the  assumptions  of  a  constant-property  fluid  and 
uniform  wall  mass  transfer  and  are  presented  for  radial 
Reynolds  numbers  (Reynolds  numbers  based  on  the  velocity  of 
the  flow  through  the  tube  wall  and  the  tube  diameter)  ranging 
from  -4.618  to  20.  Using  numerical  solutions  to  the  tube  flow 
and  energy  equations,  he  found  that  the  friction  coefficient, 
f,  was  a  function  of  both  the  radial  and  axial  Reynolds  numbers 
Re*,  and  Re«,  respectively.  For  large  values  of  Re„,  the  values 
of  the  product  f*Re  asymptotically  approached  the  value  19.739. 
He  also  found  that  the  wall  friction  decreased  with  increased 
extraction  radial  Reynolds  number  and  that  for  a  given  axial 
Reynolds  number  fluid  injection  caused  an  increase  in  the 
wall  friction. 

Yuan  and  Flnkelsteln  (11)  studied  the  effects  of  uniform 
wall  injection  and  extraction  on  two-dimensional  laminar  in¬ 
compressible  flow  in  a  tube  with  porous  walls.  They  used  a 
perturbation  method  to  solve  the  Navler-Stokes  equations  in 
cylindrical  coordinates,  and  they  observed  that  fluid  injec¬ 
tion  at  the  wall  accelerated  the  main  stream  velocity  and 
thus  increased  the  velocity  gradient  at  the  wall.  The  wall 
frictional  coefficient  was  therefore  found  to  increase  with 


injection  and  decrease  with  suction. 
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Kinney  and  Sparrow  (6)  used  an  analytical  model  to  In¬ 
vestigate  the  e££ects  o£  surface  suction  on  friction,  heat 
transfer,  and  mass  transfer  for  turbulent  flow  inside  a  tube. 
They  found  that  the  shear  force  was  greater  when  suction  was 
present  than  when  it  was  absent  and  that  high  Reynolds  number 
flows  are  more  sensitive  to  suction  than  low  Reynolds  number 
flows.  They  also  observed  that  in  the  presence  of  suction 
the  momentum  change  of  the  flow  tended  to  increase  the  pres¬ 
sure  in  the  flow  direction,  while  the  wall  shear  stress  tended 
to  decrease  the  pressure. 

For  his  doctoral  dissertation  at  the  Air  Force  Institute 
of  Technology  in  1986,  Bowman  numerically  and  experimentally 
investigated  the  effects  of  mass  injection  and  extraction  on 
the  wall  frictional  coefficients  (2).  Bowman  developed 
functional  relationships  for  friction  coefficients  to  be  used 
in  one -dimensional  heat  pipe  vapor  models  to  design  heat  pipes 
for  compressible-flow  situations  with  large  mass  injection  and 
extraction  rates.  Bowman  used  a  simulated  heat  pipe  with  air 
as  the  working  fluid  and  found  that  the  friction  coefficient 
in  general  is  a  function  of  Re%*,  Re«,  Mach  number  and  the  ratio 
of  pipe  length  to  diameter.  For  laminar  flow,  he  found  the 
relationship  to  be 

f*Re„  -  16(1.2337  -  0.2337  exp( 0 . 0363Re„) ) 

x  exp(6Ma*/5)  (1.1) 


and  for  turbulent  flow 


f  -  ft  1  ♦  55  Re°  ‘Mv/u)®  ■•(2L./D)°*  *-exp(6Ma*/5)  1  (1.2) 


where 


f*  -  0 . 046/Re*'®  (1.3) 

where  v  Is  the  radial  velocity,  u  is  the  average  axial  veloc¬ 
ity,  D  is  the  pipe  inside  diameter,  and  L.  is  the  length  of 
the  condenser . 

QtolfidlYfi.  and  Scobs 

The  main  objective  of  this  research  is  to  verify  Bowman's 
friction  coefficient  expressions  for  intermediate  radial 
Reynolds  number  compressible  flow.  To  accomplish  this,  a  sim¬ 
ulated  heat  pipe  and  a  numerical  simulation  were  used.  The 
axial  pressure  variation  was  experimentally  measured  along 
the  simulated  heat  pipe.  This  pressure  distribution  was  used 
as  input  to  a  numerical  model  designed  to  model  one-dimensional 
compressible  and  incompressible  flow  in  a  tube  with  mass  in¬ 
jection  and  extraction.  The  numerical  model  utilized  the 
pressure  distribution  to  solve  for  the  friction  coefficient 
distribution  along  the  pipe.  For  axial  flow  Hach  numbers  less 
than  0.01,  Newton's  second  law  and  incompressible  flow  equa¬ 
tions  were  used  to  solve  for  the  flow  properties;  while  for 
Mach  numbers  greater  than  0.01,  compressible  flow,  the  con¬ 
tinuity  equation  and  Shapiro's  method  of  influence  coeffi¬ 
cients  (9)  were  employed.  The  numerical  results  were  then 
compared  to  the  values  predicted  by  Bowman's  expressions. 


The  purpose  of  the  experimental  portion  of  this  study 
was  to  obtain  data  on  the  pressure  distribution  for  flow 
inside  a  simulated  heat  pipe  at  intermediate  radial  Reynolds 
numbers.  This  was  accomplished  by  means  of  a  simulated  heat 
pipe  modeled  after  the  one  used  by  Bowman  (2)  in  his  research. 
The  sample  pressure  distributions  were  then  used  as  input  to 
a  numerical  model/  to  be  discussed  in  Chapter  III/  in  order 
to  calculate  the  friction  coefficient  distribution.  This 
chapter  discusses  how  the  heat  pipe  vapor  flow  was  simulated 
and  the  experimental  configuration  used  to  model  the  flow. 

Heat  Pipe  Simulation 

A  simulated  heat  pipe  was  used  because  of  its  advantages 
over  an  actual  heat  pipe  in  experimental  studies.  The  design 
consisted  of  two  porous  pipes  of  equal  length  placed  end  to 
end  with  air  injected  through  the  walls  of  one  pipe  to  model 
the  evaporator  and  extracted  through  the  walls  of  the  other 
pipe  to  model  the  condenser.  Wageman  and  Ouevara  (10)  and 
Quale  and  Levy  (8)  showed  that  this  system  simulates  the  vapor 
dynamics  in  a  heat  pipe  quite  well.  One  advantage  of  the  sim¬ 
ulated  heat  pipe  is  that  it  is  much  simpler  to  build  and  test. 
Heat  pipes  typically  operate  at  low  pressures  with  very  small 
axial  pressure  gradients,  making  it  difficult  to  measure  pres¬ 
sures.  The  porous  pipe  system  eliminates  this  problem  since 
it  operates  at  higher  pressures  with  larger  pressure  gradients. 


yet  still  accurately  simulates  flow  in  an  actual  heat  pipe,  in 
addition,  this  system  does  not  require  complicated  heating  and 
cooling  equipment  nor  does  it  require  a  liquid  return  mecha¬ 
nism  (similar  to  a  wick  in  a  heat  pipe). 

BXPCllaentttl  Apparatus 

As  was  mentioned  earlier,  the  experimental  set  up  used 
in  this  investigation  was  modeled  after  that  used  by  Bowman. 
Much  of  the  same  equipment  he  used  was  used  for  this  experi¬ 
ment  as  well;  for  this  reason,  only  a  brief  dlscription  of 
the  set  up  will  be  given.  The  reader  is  referred  to  Bowman's 
dissertation  for  additional  details. 

In  order  to  obtain  intermediate  radial  Reynolds  numbers 
a  longer  porous  pipe  with  a  smaller  diameter  than  that  used 
by  Bowman  was  required.  A  polyethylene  porous  pipe  produced 
by  Porex  Technologies  with  a  nominal  pore  size  of  35  microns, 
a  5/8  inch  inside  diameter  and  3/4  inch  outside  diameter  was 
the  smallest  diameter  pipe  found.  Since  this  pipe  was  only 
slightly  smaller  in  diameter  than  Bowman's,  the  length  of  the 
test  section  was  increased  to  approximately  seven  feet.  The 
pipe  selected  was  only  available  in  four  foot  lengths;  there¬ 
fore,  two  four  foot  pipes  were  chosen.  One  three  inch  section 
was  cut  from  both  ends  of  each  pipe  for  calibration  tests 
leaving  two  3-1/2  foot  pipes.  One  pipe  was  used  to  model  the 
evaporator  and  was  enclosed  in  a  manifold  designed  and  built 
to  supply  the  air  for  Injection.  The  Injected  air  was  passed 
through  a  filter  just  prior  to  entering  the  test  configure- 


tion.  The  other  pipe,  used  to  model  the  condenser,  was  left 
open  to  the  atmosphere.  A  diagram  of  the  set  up  is  shown  in 
Figure  2.1.  The  two  pipes  were  placed  end  to  end  so  that  the 
air  injected  into  the  first  porous  pipe  would  flow  into  the 
second  pipe  where  it  would  be  extracted  through  the  walls  of 
the  pipe  into  the  atmosphere. 

In  the  original  test  configuration,  twelve  pressure  taps 
were  installed,  six  evenly  spaced  in  each  pipe.  When  pressure 
measurements  were  made,  it  was  found  that  the  pressure  varied 
along  the  tube  more  than  expected  in  the  condenser  section; 
so  15  additional  pressure  taps  were  added  to  that  section.  The 
final  configuration  is  shown  in  Figure  2.1.  Twenty-seven 
pressure  taps  were  used  in  the  porous  pipes,  six  in  the 
evaporator  section  spaced  six  Inches  apart  and  twenty-one  in 
the  condenser  section  spaced  every  two  inches.  In  addition, 
five  other  pressure  taps  were  used  in  the  system.  Two  were 
used  to  measure  the  pressure  around  the  orifice  plate,  one 
above  the  plate  and  one  below  the  plate.  One  tap  was  used 
to  measure  the  pressure  at  the  upstream  end  of  the  air  supply 
manifold  and  two  taps  were  used  to  measure  the  pressure  at 
the  downstream  end  of  the  manifold.  Stainless  steel  0.060 
inch  outside  diameter  tubes  were  used  as  pressure  taps.  The 
tubes  were  Inserted  and  cemented  into  0.058  inch  diameter 
holes  in  the  porous  pipe  to  produce  a  tight  fit.  The 
stainless  steel  tubes  were  Installed  such  that  their  ends 
were  flush  with  the  inside  wall  of  the  porous  pipe. 


Figure  2.1  Experimental  Test  Configuration 
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The  calibration  of  the  poroua  pipe  was  made  using  a 
relationship  developed  by  M.  Muskat  (7:66).  Muskat  found  that 
the  mass  flux  through  a  porous  medium  is  related  to  the  pres¬ 
sure  difference  across  the  porous  medium  by  the  following 
expression: 


(Apa)  *  A(pv ) *  +  B(pv) 


(2.1) 


Bach  of  the  four  porous  pipe  samples  taken  from  the  pipe 
used  in  the  test  configuration  were  installed  into  the  test 
manifold  designed  by  Bowman  (Figure  2.2).  For  each  sample, 
APa  and  pv  were  measured  at  various  mass  flow  rates.  The 
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Side  View 


Top  View 


Pigure  2.2  Porous  Pipe  Calibration  Manifold 

(Bowman  2: 29) 


curve  was  expected  to  be  very  close  to  linear  at  very  low 
nass  flow  rates,  with  the  contribution  of  the  A(pv)*  tern 
being  quite  snail  compared  to  the  B(pv)  tern.  This  linear 
behavior  was  observed  in  that  region,  however  the  curve  was  off¬ 
set  from  the  origin  by  the  amount  AP*  *  126,921.84  lbfa/ft4. 
These  results  are  shown  in  Pigure  2.3.  This  offset  in  the 
curve  was  probably  due  to  the  tendency  of  the  pressure  trans¬ 
ducer  in  the  Scanivalve  to  read  the  pressures  higher  than 
what  they  actually  were.  The  atmospheric  pressure  readings 
Mde  by  the  Scanivalve  were  consistently  higher  than  those 
taken  from  the  barometer  in  the  laboratory.  To  correct  for 
this  effect,  the  offset  value  was  subtracted  from  each  AP* 
reading.  The  values  for  the  constants  A  and  B  in  Equation 
(2.1)  were  calculated  using  a  least  squares  technique  to  curve 
fit  the  corrected  calibration  data  and  were  found  to  be 
4.679  x  10*  ftVsec*  and  4.2843  x  10T  lbf/ft-sec,  respectively. 
The  final  calibration  results  are  shown  in  Pigures  2.4  and 
2.5.  Pigure  2.5  is  an  enlargement  of  the  region  of  the  graph 
in  Figure  2.4  for  the  low  mass  flow  rates.  Equation  (2.1) 
was  needed  for  the  numerical  model  to  determine  the  mss  flux 
for  a  given  pressure  differential  across  the  pipe  walls. 

QAtA.  BfiflgmeaftDt  and  Acquisition 

A  Zenith  Z-100  computer  was  used  for  the  data  acquisition. 
Pressure  MasureMnts  were  Mde  using  a  T-type  Scanivalve, 
from  Scanivalve  Corporation,  with  a  model  PDCR42  differential 
pressure  transducer  having  a  pressure  range  of  -100  to  100 


Porous  Pipe  Calibration  Data  Showing  Offset  From  Origin 


Figure  2.4  Porous  Pipe  Calibration  Results 


2.5  Enlargement  of  A(P  )  Range  0  -  500,000  lbf  /ft  in  Calibration  Results 


psig.  A  Dual  Systems  control  model  AIH12  analog-to-dlgltal 
card  In  the  computer  was  used  to  digitize  the  pressure  trans¬ 
ducer  Input  into  the  computer.  The  offset  option  on  the  card 
and  amplification  were  set  to  read  signals  from  -100  to  +100 
mV  with  a  resolution  of  0.048  mV.  Fifty  data  samples  were 
averaged  for  each  point  to  reduce  noise.  For  mass  fluxes  above 
0.035  lbf-sec/ft*,  a  square-edged  orifice  plate  was  used  to 
measure  the  total  mass  flow  rate  entering  the  test  configura¬ 
tion.  The  pressure  differential  across  the  orifice  plate  was 
measured  using  both  the  Scanivalve  and  a  water  manometer.  The 
mass  flow  rate  was  calculated  using  the  following  expression 
(1): 


m  -  0. 52502  (CYd*F. )  <  p*  (p*  -  p»)x'V(l  -  P*) 


(2.2) 


where 


C  -  K/B  -  K(  1  -R,  )*'*  (2.3) 

K  *  K.(  1  ♦  A/Re*)  (2.4) 

K.  «  Kk. (  10*d/(10«d  ♦  15A)  )  (2.5) 

Ka.  -  .5993  «■  . 007/D  ♦  (.364  ♦  .076/Dx'*)P* 

♦  .4(1.63  -  1/D) ( ( .07  *■  .5/D)  -  P*)*'» 

-  (.009  ♦  .034/D) (.5  -  P*)*'* 

♦  (65/D*  +  3)<p*  -  .7)*'*  (2.6) 

A  -  d( 830  -  5000P  ♦  9000p*  -  4200P*  ♦  530/Dx'*)  (2.7) 

Y  •  1  -  (.41  ♦  .  35P*4  )  (X/k )  (2.8) 

X  -  P/Pi  (2.9) 


Re*  -  48m/  n  p  d 


(2.10) 


Very  low  mass  flow  rates,  such  as  those  measured  In  the  cali¬ 
bration  tests,  were  measured  using  a  flowmeter  with  a  range 
of  0  to  0.44  cfm.  The  pressure  was  measured  at  a  point  just 
upstream  of  the  flowmeter  in  order  to  calculate  the  density 
of  the  air.  In  addition,  a  valve  was  installed  just  down¬ 
stream  of  the  flowmeter  to  vary  the  upstream  pressure  In  order 
to  read  higher  mass  flow  rates  with  the  meter.  It  was  impor¬ 
tant  to  get  accurate  readings  in  the  low  mass  flow  rate  region 
because  there  were  several  sections  of  the  simulated  heat  pipe 
where  the  lnjectlon/suction  rates  were  quite  small.  Bquatlon 
(2.1)  with  the  constants  A  and  B  calculated  from  data  for 
mss  fluxes  only  above  0.04  lbf-sec/ft*  under  predicted  the 
mass  fluxes  at  low  Apa« 

Experimental  Results 

The  pressure  variation  obtained  along  the  pipe  was  a 
result  of  two  effects:  the  friction  encountered  along  the 
pipe  walls  and  the  rate  of  change  in  momentum  along  the  pipe. 
The  total  effect  of  the  friction  force  was  manifested  in  the 
pressure  drop  between  the  two  ends,  or  the  lack  of  complete 
pressure  recovery.  The  pressure  drop  due  to  friction  was 
found  to  Increase  with  increased  mass  flow  rate.  Figure  2.6 
shows  the  pressure  distribution  for  five  different  mass  flow 
rates.  Additional  data  for  other  mass  flow  rates  can  be 
found  in  Appendix  A.  As  mass  was  added  in  the  evaporator 
region,  0  <  X/L  <  0.5,  the  fluid  accelerated  with  increasing 
X/L  with  a  corresponding  drop  in  pressure.  The  mass  removal 


in  the  condenser  section  caused  the  £lov  to  decelerate  and 
the  pressure  to  rise. 
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This  chapter  describee  the  numerical  Model  used  to  sim¬ 


ulate  the  flow  generated  by  the  experimental  model.  The 
governing  equations  used  to  solve  for  the  flow  properties  will 
be  discussed  first,  followed  by  a  dlscription  of  the  solution 
technique  used.  Finally,  the  results  obtained  will  be  pre¬ 
sented  . 

Governing  Equations. 

All  of  the  flowfields  generated  with  the  experimental 
configuration  contained  sizable  regions  of  very  low  Hach  num¬ 
ber  (Ma  <  0.1).  Therefore,  two  models  were  used  in  the  num¬ 
erical  simulation:  an  incompressible  model  for  Hach  numbers 
less  than  0.01,  and  a  compressible  model  for  Hach  numbers 
greater  than  0.01.  One-dimensional,  steady,  adiabatic  flow 
was  assumed  for  both  models. 

Incompressible  Hodel .  The  incompressible  model  consisted 
primarily  of  Newton's  second  law  and  incompressible  flow  rela¬ 
tions.  To  approximate  the  shear  stress  at  the  pipe  wall, 
Newton's  second  law  was  applied  to  an  element  of  fluid  in  the 
pipe  (Figure  3.1). 


Figure  3.1  Blement  of  Fluid  for  Equation  of  Hotion 


Applying  Newton's  second  law  for  steady  flow  to  the  fluid 
element  gives 

•aUl/g.  -  »iUx/g-  ■  (pxAx  -  Pa  Aa )  -  TWA»  (3.1) 

where  A»  is  the  wall  surface  area  for  the  increnent  ax. 

Equation  (3.2)  can  be  solved  for  Tw  giving 

T «r  »  (  ItiUi/g.  -  MiUa/g.  ♦  (pa.  -  Pa) A-  ]/A*  (3.2) 

where  A.  is  the  cross-sectional  area  equal  to  Ax  and  Aa.  The 
expression  for  the  friction  coefficient  is  given  by 

f  “  2  T*,  g./  pu*  (3.3) 

where 

u  -  (ux  ♦  u.)/2  (3.4) 

Substituting  the  expression  for  t«#  given  by  Equation  (3.2)  in 
Equation  (3.3)  gives  an  expression  for  the  friction  coefficient 
in  terse  of  quantities  which  can  be  calculated  fro*  the  input 
pressure  distribution 

f  ■  2t  AM  ♦  (px  -  pa) A.  lg./  puaA»  (3.5) 

where 

•  • 

AM  ■  ixUx/g.  -  »aUa/g.  (3.6) 

The  axial  velocities  along  the  pipe  were  calculated  fro* 
a  ease  balance  on  an  element  of  fluid  in  the  pipe  (Figure  3.2). 
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Figure  3.2  Element  of  Fluid  Used  For  Hass  Balance 


Performing  the  mass  balance  on  the  fluid  element  in  Figure  3.2 


PlUaAa  ■  PlUlAs  —  (pV)awAy 


(3.7) 


where  (pv)av«  is  the  average  mass  flux  through  the  pipe  walls 
over  the  Increment  ax  found  by  applying  Equation  (2.1)  at 
stations  1  and  2  (shown  in  Figure  3.2)  and  taking  the  average 
of  the  two.  The  negative  sign  on  the  quantity  (pv)av«  is  a 
result  of  the  sign  convention  on  the  mass  flux.  For  injection, 
(pv)  is  considered  negative  and  considered  positive  for  suction. 
Equation  3.7  can  be  solved  for  Ua  giving 


Ua  •  piUi/Pa  “  (pv)avaAf/Ptka 


(3.8) 


The  velocity  ux  is  taken  from  the  previous  step  and  its  ini¬ 
tial  value  is  zero. 

The  pressure  at  each  step  was  found  by  interpolating 
from  the  input  pressure  distribution  using  a  second  degree 
interpolation  formula.  Figure  3.3  illustrates  this  interpo¬ 
lation  procedure. 


M 


Figure  3.3  Pressure  Interpolation 
Refer ing  to  Figure  3.3  for  xi  <  x  <  x*,  p(x)  Is  given  by 

P(x)  -  pxK*  ♦  p*Ka  ♦  p>K>  (3.9) 

vhere 

Kx  -  ( (X  -  Xa ) ( X  -  X*)l/l(xx  “  X*)(Xx  -  Xa ) ) 

K»  •  ((x  -  XxHx  -  Xa )  )/(  (Xa  -  Xx)(X*  -  X|)1  (3.10) 

K»  ■  1(X  -  XxXx  -  Xa))/[(Xj  -  Xx)(Xa  -  Xa)] 

The  temperature  (T),  viscosity  (p),  axial  Reynolds  num¬ 
ber  (Re-),  and  radial  Reynolds  number  (Rev)  were  calculated 


from  the  following  expressions: 

T  -  T./[  1  ♦  (k  -  1 )Ma*/2  ]  (3.11) 

M  -  2.271-08  (T)  ■/’*/  (T  ♦  198.6)  (3.12) 

Re-  ■  4  m.vw  /(n  d  p  g.)  (3.13) 

Re-  -  (pv )— D  / (  p  g- )  (3.14) 


where  T  in  Equation  (3.11)  is  in  degrees  Fahrenheit,  and 


Is  the  average  of  mx  and  m*.  Equations  (3.5),  (3.6)  -  (3.14), 
and  (2.1)  thus  comprise  the  equations  for  the  Incompressible 
model . 

rnfr>«gible  Model .  The  compressible  model  consisted  of 
Equations  (3.9)  -  (3.14)  along  with  the  differential  form  of 
the  continuity  equation  and  Shapiro's  method  of  Influence 
coefficients  for  Mach  number  variation  with  friction  and  wall 
mass  transfer  to  solve  for  the  velocity  and  friction  coeffi¬ 
cient,  respectively.  The  differential  form  of  the  continuity 
equation  in  terms  of  pressure,  temperature  and  Mach  number  is 
given  by 

da  -  da  +  dMfl.  +  1  dX  (3.15) 

m  ”  P  Ma  2  T 

Using  the  equation  of  state  for  an  Ideal  gas,  p  *  p/RT,  and 
Equation  (3.11)  in  equation  3.15  gives  the  differential  form  of 
the  continuity  equation  in  terms  of  pressure  and  Mach  number. 

d(Mal»  (3.16) 

Ma* 

Equation  (3.16)  can  then  be  integrated  to  solve  for  Maa  ,  which 
is  given  by  the  following  expression 

Ma»2  -  (-1  ♦  11  ♦  2 ( k  -  l)C)>/(k  -  1)  (3.17) 

where 

C  -  Max2  (1  ♦  (k  -  1  )MaxZ/2  l(m*/mx) */(p*/px)*  (3.18) 


da  =  da  +  1 

m  p  2 


1  ±  (K  - 

1  ♦  (k  -  l)Ma*/2 

m 
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The  equation  for  Mach  nuaber  variation  with  friction 
and  wall  aass  transfer  utilizing  Shapiro's  method  of  influence 
coefficients  is  given  by 

•  • 

d(Ma* )/Ma*  -  F*4  £  dx/D  +  P«  da/a  (3.19) 

where  F«  and  P«  are  the  influence  coefficients  defined  by 

Vm  -  Ma*(  1  ♦  (k  -  l)Ma*/2 1/(1  -  Ma*)  (3.20) 

Pw  «  2(1  ♦  kMa* ) ( 1  +  (k  -  l)Maa/2]/(l  -  Maa)  (3.21) 

where  x  is  the  axial  coordinate,  equation  (3.19)  can  be  rear¬ 
ranged  into  the  following  fora 

(1  -  Ma* )  dHaa  4fdx  ,  2(1  +  kMaa)  da  (3.22) 

kMa*( 1  ♦  (k  -  1 )Ma*/2 1  D  kMa*  a 

Equation  (3.22)  is  in  a  fora  which  can  be  Integrated  to  yield 
an  expression  for  the  friction  coefficient.  After  perforalng 
the  integration  and  rearranging,  the  following  expression 
results 

f  «  DC  In (Max2 /Mas2  )  -  C-  ln(p./px ) J/( 2ax)  (3.23) 

where 

C«  -  (1  ♦  kMa* )/( kMa* )  (3.24) 

ite*  -  (Max2  ♦  Ma*z)/2  (3.2$) 

and  where  Max2  Is  found  froa  Equation  (3.17)  and  Max2  is  taken 
froa  the  previous  step.  Equations  (3.9)  -  (3.15),  (3.17), 
(3.23)  and  (2.1)  fora  the  set  of  equations  used  in  the  coapres 


slble  aodel. 


Solution  Method 

A  marching  technique  was  used  to  solve  for  flow  proper¬ 
ties  along  the  simulated  heat  pipe.  The  pressure  lnforstttlon 
obtained  from  the  experimental  work  was  first  read  in  to  the 
program;  then  the  various  constants  for  the  pipe  geometry  and 
the  air  properties  and  the  variables  at  the  upstream  end  of 
the  pipe  were  set.  The  set  of  equations  for  the  incompres¬ 
sible  model  were  then  solved  to  find  the  Kach  number,  axial 
Reynolds  number,  and  friction  coefficient  for  each  ax  incre¬ 
ment  along  the  pipe  until  a  Mach  number  of  0.01  was  reached. 
After  that,  the  compressible  model  equations  were  used  to  de¬ 
termine  the  properties.  Sample  calculations  for  the  first 
iteration  sade  by  the  program  are  given  in  Appendix  B  and  a 
copy  of  the  code  used  is  included  as  Appendix  C. 

To  verify  the  accuracy  of  the  numerical  simulation,  three 
checks  were  made  in  the  program.  The  first  check  was  a  com¬ 
parison  of  the  mass  flow  rate  entering  the  condenser  with  the 
mass  flow  rate  entering  the  air  supply  manifold  for  the  evap¬ 
orator  to  see  if  the  two  were  equal.  The  second  check  was  a 
force  balance  over  the  entire  pipe.  When  considering  the  pipe 
as  a  whole,  the  only  forces  acting  on  the  fluid  are  the  pres¬ 
sures  at  the  two  ends  and  the  total  wall  shear  force  (see 
Figure  3.4). 


v*> 


Figure  3.4  Force  Balance  on  Entire  Heat  Pipe 

A  summation  of  these  forces  gives 

IF  -  (p»  -  p-)A.  -  T«#  Ag»  (3.26) 

where  A»  is  the  surface  area  cf  an  increment  ax  of  pipe. 
Ideally,  this  summation  should  be  zero,  since  the  drop  in 
pressure  between  the  two  ends  of  the  pipe  is  a  result  of  the 
wall  shear  stress.  The  last  check  made  was  a  mass  balance  over 
the  entire  pipe  to  see  if  the  total  mass  injected  into  the  pipe 
equaled  the  total  mass  removed  from  the  pipe. 

Muaeilcai  ReaulU 

All  of  the  numerical  simulations  made  using  the  code  in 
Appendix  C  modeled  the  porous  pipe  system  used  in  the  exper¬ 
imental  study  for  various  supply  mass  flow  rates.  The  results 
are  plotted  as  the  friction  coefficient  -  axial  Reynolds  num¬ 
ber  product  (f*Re)  versus  Mach  number  for  the  evaporator  in 
Figure  3.S  and  for  the  condenser  in  Figure  3.6. 


Effect  of  Compressibil 


Figure  3.6  Effect  of  Compressibility  in  Condenser 


IV.  Discussion  q±  Results 


This  chapter  discusses  the  reaulta  obtained  in  both  the 
experimental  and  numerical  portions  of  this  research.  The 
reader  is  referred  to  Figure  2.6  for  the  discussion  of  the 
experimental  results  and  to  Figures  3.5  and  3.6  for  the  dis¬ 
cussion  of  the  numerical  results. 

BXPOElBCntal  Results  Discussion 

An  important  part  of  the  experimental  study  was  the 
calibration  of  the  porous  pipe.  It  was  particularly  impor¬ 
tant  to  obtain  accurate  data  for  the  very  low  mass  flux  range, 
0.0  -  0.005  lbf-sec/f t*,  since  the  simulated  heat  pipe  oper¬ 
ated  with  several  regions  of  very  low  mass  flux.  The  four 
pipe  samples  taken  from  the  porous  pipe  used  in  the  test  con¬ 
figuration  were  originally  calibrated  using  the  orifice  plate- 
water  manometer  system  (described  in  Chapter  II)  to  measure 
the  mass  flow  rate  entering  the  test  sample.  This  system  only 
allowed  measureamnts  of  mass  flow  above  approximately  0.03 
lbf-sec/ft*,  since  the  ratio  of  the  orifice  plate  bore  to  the 
air  supply  pipe  diameter  was  0.327  and  values  of  the  discharge 
coefficient  are  not  available  for  ratios  below  0.3.  Curve- 
fitting  the  data  from  these  first  calibration  tests  produced 
values  for  the  constants  A  and  B  in  Equation  (2.1)  which  under 
predicted  the  mass  flux  for  small  values  of  Ap*.  Further  cal¬ 
ibration  tests  were  then  made  on  two  of  the  four  samples,  one 
sample  from  the  pipe  used  in  the  evaporator  section  and  one 


froa  the  pipe  used  in  the  condenser  section,  using  a  flows* ter 
with  a  range  of  0  -  0.44  cfa.  Only  two  of  the  samples  were 
calibrated  because  the  other  two  were  too  dirty  and  clogged 
froa  the  first  set  of  calibration  tests  to  yield  accurate  re¬ 
sults.  The  data  obtained  froa  the  second  set  of  calibration 
tests  exhibited  the  linear  behavior  expected  in  the  low  aass 
flux  region,  but  the  curve  showed  an  offset  froa  the  origin, 
where  it  should  have  passed  through.  This  offset  probably 
resulted  froa  a  combination  of  two  things.  One  is  that  the 
pressure  transducer  tended  to  give  pressure  readings  that  were 
too  high.  This  was  observed  when  the  transducer  readings  for 
the  ataospherlc  pressure  were  compared  to  the  readings  taken 
froa  the  barometer  in  the  laboratory,  which  is  considered  the 
more  accurate  of  the  two.  The  transducer  consistently  gave 
higher  readings  than  the  barometer,  by  approximately  IS.  The 
other  possible  reason  for  the  offset  in  the  curve  is  the  fact 
that  the  ataospherlc  pressure  used  to  calculate  the  Ap*  values 
was  taken  froa  the  barometer,  not  froa  the  pressure  transducer. 
Therefore  in  each  of  the  Ap*  teras,  the  pressure  value  for 
the  inside  of  the  pipe  was  too  high  and  the  external  pressure 
value  was  fairly  accurate  thus  making  the  Ap*  values  higher 
than  they  actually  were.  As  was  stated  in  Chapter  II,  the 
offset  was  corrected  for  by  subtracting  the  amount  of  the 
offset  from  each  of  the  Ap*  readings.  The  constants  obtained 
using  the  additional  calibration  data  yielded  values  for  the 
constants  A  and  B  that  improved  the  results  in  the  low  mass 
flux  region.  Upon  examining  the  range  of  mass  fluxes  in  the 


simulated  heat  pipe,  it  was  found  that  most  of  the  mss  flux 
values  were  below  0.05  lbf-sec/ft*.  Another  curve-fit  was 
then  wide  of  only  the  calibration  data  up  to  0.05  lbf-sec/ft», 
which  produced  values  for  A  and  B  that  further  improved  the 
results  in  the  low  mss  flux  region.  Since  the  Mils  of  the 
pipe  used  in  these  experiments  were  approxiMtely  half  as 
thick  as  those  used  in  Bowman's  experiments,  the  value  of  B 
was  expected  to  be  less  than  the  value  that  Bovnan  calculated. 
As  it  turned  out,  this  value  was  about  one  fourth  that  of 
Bowun '  s . 

The  pressure  distributions  obtained  from  the  tests  with 
the  simulated  heat  pipe  showed  a  fluctuation  in  the  pressure 
distribution  in  the  condenser  section  rather  than  a  smooth  in¬ 
crease  of  pressure  from  the  beginning  to  the  end  of  the  conden¬ 
ser.  Kinney  and  Sparrow  (6)  noted  that  in  the  presence  of  wall 
suction,  the  change  in  momentum  of  the  flow  tended  to  increase 
the  pressure  in  the  flow  direction,  while  the  wall  shear 
stress  tended  to  decrease  the  pressure.  In  the  first  half 
of  the  condenser,  the  change  in  momentum  was  very  bmII  due 
to  the  smII  amount  of  suction  through  the  Mils  and  ms  ap¬ 
proximately  the  same  order  of  Mgnitude  as  the  Mil  shear 
stress.  Therefore  the  momentum  change  was  not  the  dominant 
effect  in  this  region  and  consequently  the  pressure  did  not 
exhibit  a  smooth  increase. 

Due  to  the  compressible  nature  of  the  flow  and  the  low 
design  pressure  drop  across  the  pipe  Mils,  the  mss  injec¬ 
tion  and  removal  along  the  pipe  studied  was  non-uniform.  As 


was  discussed  previously,  the  Mass  flux  through  the  pipe  wall 
is  dependent  on  the  pressure  difference  across  the  pipe  wall. 
Since  the  external  pressure  was  relatively  constant  and  the 
internal  pressure  varied  axially  along  the  pips,  the  mss 
injection  and  removal  also  varied.  The  thinness  of  the  pipe 
walls  Mde  the  pressure  drop  across  the  pipe  low,  thus  pre¬ 
senting  less  resistance  to  the  mss  flux.  The  pipe's  inter¬ 
nal  pressure  consequently  reulned  quite  close  to  the  external 
pressure  for  a  large  portion  of  both  the  condenser  and  the 
evaporator  sections.  The  pressure  dropped  off  sharply  in  the 
last  six  inches  of  the  evaporator,  where  Most  of  the  mss  ad¬ 
dition  occurred,  and  rose  sharply  in  the  last  half  of  the  con¬ 
denser,  where  Most  of  the  mss  was  reMoved.  This  suggests 
that  there  My  be  an  optiMUM  length  or  L/D  for  the  pipe. 

Muaeilcal  Rcauita  Discussion 

In  all  of  the  nuMrical  simulations,  the  flow  in  the 
evaporator  section  was  assumed  to  be  laminar  throughout.  This 
was  a  reasonable  assumption  since  Bowsan  showed  that  the  flow 
reMins  laminar,  when  injection  is  present,  up  to  very  large 
axial  Reynolds  numbers  (on  the  order  of  10*).  The  values  of 
the  friction  coefficients  obtained  from  the  nuMrical  simula¬ 
tions  were  somewhat  higher  than  the  values  predicted  by 
BowMn's  expression.  In  the  evaporator  section,  where  mss 
Injection  was  present,  Kinney's  (5)  results  Indicate  that  for 
laminar,  incompressible  flow  with  large  injection  wall  Reynolds 
numbers,  the  product  f*Ren  should  approach  approxiMtely  20. 


Por  the  incompressible  region  of  the  evaporator  used  In  this 
study,  the  calculated  friction  coefficients  are  approximately 
one  order  of  magnitude  greater  than  Kinney's  results.  Bowman's 
laminar  flow  friction  coefficient  expression  included  the  ef¬ 
fect  of  the  compressibility  of  the  flow  and  thus  extended 
Kinney's  results  for  higher  Mach  numbers.  Again  the  calculated 
friction  coefficients  are  approximately  one  order  of  magnitude 
greater  than  the  values  predicted  by  Bowman's  expression. 
Kinney's  and  Bowman's  results  indicate  that  the  product  f*Re 
remains  essentially  constant  (approximately  20)  for  low  Mach 
numbers  and  Bowman's  expressions  predict  an  Increase  as  the 
Mach  number  increases  toward  1  and  greater.  The  results  from 
this  research  ware  found  to  exhibit  this  trend  for  the  f*Rea 
product  versus  Mach  number,  which  can  be  seen  in  Pigure  3.5, 
only  the  Increase  with  Mach  number  begins  at  lower  Mach  num¬ 
bers  than  Bowman's  results. 

The  flow  in  the  condenser  section  was  assumed  to  be  tur¬ 
bulent  throughout  in  all  of  the  numerical  simulations.  Bowman 
showed  that  once  the  flow  entered  the  condenser,  where  an  ad¬ 
verse  pressure  gradient  was  present,  it  transitioned  to  tur¬ 
bulent.  He  also  found  that  the  transition  always  started  at 
the  beginning  of  the  condenser  near  the  pipe  wall,  when  the 
flow  was  subsonic,  and  then  spread  toward  the  center  of  the 
pipe  as  the  flow  progressed  through  the  condenser.  The  loca¬ 
tion  in  the  condenser  where  the  turbulence  reached  the  pipe 
center  moved  downstream  as  the  axial  velocity  increased. 

8ince  the  transition  regions  for  the  axial  velocities  attained 


at  the  entrance  to  the  condenser  In  this  research  were  quite 
s stall  relative  to  the  overall  condenser  length,  neglecting 
the  transition  region  and  assuming  turbulent  flow  throughout 
was  a  reasonable  assumption.  The  values  computed  for  the 
friction  coefficients  in  the  condenser  section  were  lower,  by 
approximately  one  order  of  magnitude,  than  the  values  predict¬ 
ed  by  Bowman's  turbulent  flow  friction  coefficient  expression. 
Bowman's  expression  indicated  an  increase  in  the  f*Re»  product 
with  increased  Mach  number.  When  the  results  from  the  numer¬ 
ical  simulation  of  the  condenser  were  plotted  as  f*ReM  versus 
Ha,  the  trend  predicted  by  Bowman's  expression  was  not  seen 
as  prominently  as  the  trend  shown  in  the  evaporator,  but  it 
was  present  (see  Figure  3.6  ). 

There  were  several  sources  of  error  in  the  experimental 
study  which  could  account  for  the  difference  between  the 
calculated  and  expected  results.  As  was  mentioned  earlier, 
there  was  a  problem  with  the  pressure  transducer  in  the 
Scanivalve  tending  to  read  the  pressures  higher  than  what 
they  actually  were.  This  problem  primarily  affected  the  cal¬ 
ibration  results  by  giving  a  higher  pressure  differential  for 
a  given  mass  flux,  which  thus  Introduced  error  into  the  com¬ 
puted  values  of  the  constants  A  and  B  in  Bquation  (2.1).  Also 
there  were  no  pressure  taps  in  the  last  six  inches  of  the 
evaporator  section,  preceding  the  condenser,  where  the  pressure 
showed  a  sharper  drop.  Therefore,  the  accuracy  of  the  inter¬ 
polated  pressures  in  that  region  is  questionable.  In  all  the 
numerical  simulations,  the  friction  coefficients  near  the  end 


of  the  evaporator  were  negative  indicating  that  the  change  in 
momentum  across  the  Increment  was  greater  than  the  change  in 
pressure.  This  meant  that  the  calibration  equation  (Equation 
(2.1))  Was  giving  mass  flux  values  that  were  too  high  for  the 
given  pressure  difference/  thus  indicating  that  the  calibra¬ 
tion  constants,  A  and  B,  were  inaccurate. 

The  checks  made  in  the  numerical  simulation  code  to  ver¬ 
ify  the  accuracy  of  the  calculations  all  indicated  that  the 
siaulatlons  were  giving  accurate  results.  The  force  balance 
check  yielded  quite  small  positive  values  in  all  of  the  simu¬ 
lations  and  the  mass  balance  check  yielded  small  negative 
values  for  the  low  supply  mass  flow  rates  and  small  positive 
values  for  the  higher  supply  mass  flow  rates.  In  all  of  the 
simulations,  the  mass  flow  rate  entering  the  condenser  was 
found  to  be  too  high.  To  correct  this,  the  equation  used  to 
calculate  the  mass  flux  through  the  wall  (Equation  (2.1))  was 
multiplied  by  a  correction  factor.  This  correction  factor 
Increased  with  higher  supply  mass  flow  rates  and  ranged  from 
approximately  0.70  to  0.65.  The  addition  of  the  correction 
factor  also  improved  the  force  and  mass  balances. 


This  chapter  discusses  the  conclusions  made  fro*  the 
findings  of  this  research.  In  addition,  a  brief  list  of  rec¬ 


ommendations  for  further  research  is  presented. 

Conclusions 

The  data  obtained  from  this  research  exhibited  the  sane 
basic  trends  predicted  by  previous  research,  however  differing 
by  approximately  one  order  of  magnitude.  The  change  in  momen¬ 
tum  of  the  flow  has  much  more  Influence  on  the  pressure  change 
along  the  pipe  than  the  friction  force.  Therefore,  even 
though  the  values  of  the  friction  coefficients  calculated  in 
this  work  were  much  greater  than  those  predicted  by  Bowman's 
expressions,  the  actual  difference  in  the  pressure  distribu¬ 
tions  resulting  from  Bowman's  friction  values  and  those  meas¬ 
ured  in  the  experimental  portion  of  this  research  may  not  be 
as  significant.  Another  numerical  simulation  utilizing 
Bowman's  friction  coefficient  expressions  to  solve  for  the 
axial  pressure  distribution  would  be  needed  in  order  to  make 
the  comparison.  The  friction  coefficient  calculations  were 
also  found  to  be  quite  sensitive  to  the  values  computed  for 
the  calibration  constants,  A  and  B.  The  difference  in  the 
friction  coefficients  calculated  from  this  research  and  that 
of  Bowman's  dropped  from  two  orders  of  magnitude  to  one  order 
of  magnitude  between  the  second  and  final  curve-fits  of  the 
calibration  data.  The  analysis  of  the  results  from  the  nu- 


merical  simulations  seams  to  Indicate  that  the  either  the  data 


obtained  in  the  calibration  tests  is  inaccurate  or  the  equa¬ 
tion  used  to  curve-fit  the  data  (Equation  (2.1))  does  not  qive 
an  accurate  enough  representation  of  the  calibration  results. 
Whichever  is  the  case,  the  inaccuracy  of  the  calibration  equa¬ 
tion  seems  to  be  the  primary  source  of  the  discrepency  with 
the  results  from  previous  research. 

BccomrifHiatlQjMi 

Upon  examining  the  procedures  used  in  this  study,  several 
recommendations  can  be  made  for  improvements  to  this  experi¬ 
ment  and  possibilities  for  future  research.  One  recommenda¬ 
tion  which  would  improve  the  readings  obtained  from  the  pres¬ 
sure  transducer  is  to  subtract  the  transducer's  atmospheric 
pressure  reading  from  each  reading  made  by  the  transducer  to 
get  a  gauge  pressure.  The  atmospheric  pressure  taken  from  a 
barometer  can  then  be  added  to  these  gauge  pressures  to  get 
the  absolute  pressure.  This  would  help  eliminate  the  offset 
in  the  calibration  curve,  which  was  observed  in  this  study. 
Some  of  the  recommendations  for  future  research  include: 

1)  Perform  experiments  to  measure  the  velocity  gradi¬ 
ents  along  the  wall  as  well  as  the  pressure  distribution  along 
the  pipe.  This  would  provide  more  accurate  information  for 
computing  the  shear  stress  at  the  wall  and  should  give  more 
accurate  results  for  the  friction  coefficients. 

2)  Perform  a  one -dimensional  numerical  simulation  which 
utilizes  Bowman's  friction  coefficient  expressions  to  solve 
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for  the  pressure  distribution  along  e  heat  pipe.  Then  perform 
experiments  to  measure  the  pressure  distribution  in  a  simu¬ 
lated  heat  pipe  operating  under  the  same  conditions  to  check 
the  computed  results. 


Appendix  a:  bam 


Data 


This  appendix  contains  the  raw  data  obtained  in  the 
experimental  portion  if  this  research.  Bach  file  gives  the 
mass  flow  rate  entering  the  test  section  and  34  pressure 
measurements.  Figure  2.1  shows  the  location  of  the  pressure 
taps  in  the  test  section.  In  general,  pressure  measurements 
1  and  2  were  for  the  pressure  above  and  below  the  orifice 
plate  respectively,  3  for  the  pressure  at  the  upstream  end  of 
the  evaporator  manifold,  4  and  5  for  the  downstream  end  of  the 
evaporator  manifold,  6  and  7  for  the  atmospheric  pressure, 
8-13  for  inside  the  evaporator  section  of  the  porous  pipe, 
and  14  -  34  for  Inside  the  condenser  end  of  the  porous  pipe. 
Pressure  taps  8-13  were  approximately  six  Inches  apart,  13 
and  14  eight  inches  apart,  14  -  34  two  inches  apart.  Pressure 
tap  34  was  located  at  the  end  of  the  condenser  section  of  the 
porous  pipe.  The  total  length  of  the  test  section  was  6.976 
feet.  The  air  supply  temperature  for  all  runs  was  72  °F  + 


Data  Pile:  D2209B 


Supply  Mass  Plow  Rate  (lbn/sec)  *  0.043651 


PORT 

PRESSURE (ps la) 

PORT 

PRB8SURB(psia ) 

1 

114.4222 

18 

14.4009 

2 

114.3897 

19 

14.4229 

15.9113 

20 

14.4306 

4 

15.8919 

21 

14.5212 

5 

15.8867 

22 

14.5186 

6 

14.2250 

23 

14.5302 

7 

14.2186 

24 

14.5341 

8 

15.7949 

25 

14.2302 

9 

15.7664 

26 

14.7397 

10 

15.7238 

27 

14.7526 

11 

15.6423 

28 

14.8147 

12 

15.4638 

29 

14.8406 

13 

15.1871 

30 

14.8949 

14 

14.4009 

31 

14.9220 

15 

14.3091 

32 

14.9388 

16 

14.3298 

33 

14.9647 

17 

14.3337 

34 

14.9634 

Data 

Pile:  D2209C 

Supply  Mass  Plow  Rate 

( lbn/sec )  * 

0.058069 

PORT 

PRESSURE ( ps la ) 

PORT 

PRESSURE ( ps la ) 

1 

114.3252 

18 

14.3970 

2 

114.2675 

19 

14.4229 

3 

16.9626 

20 

14.4113 

4 

16.9264 

21 

14.5431 

5 

16.9173 

22 

14.5108 

6 

14.2121 

23 

14.5509 

7 

14.2250 

24 

14.5121 

8 

16.8527 

25 

14.2289 

9 

16.8152 

26 

14.9117 

10 

16.7583 

27 

14.9634 

11 

16.6535 

28 

15.0824 

12 

16.4195 

29 

15.1302 

13 

15.9669 

30 

15.2492 

14 

14.4203 

31 

15.3048 

15 

14.3052 

32 

15.3643 

16 

14.2929 

33 

15.3733 

17 

14.2949 

34 

15.3772 

Data  Pile:  D2209D 


Supply  Hass  Plow  Rate  (lba/sec)  -  0.070994 


PORT 

PRESSURE (pa  la) 

PORT 

PRESSURE (pa la) 

1 

114.1545 

18 

14.3983 

2 

114.0680 

19 

14.4113 

3 

18.1315 

20 

14.3673 

4 

18.0876 

21 

14.5522 

5 

18.0501 

22 

14.4746 

6 

14.2186 

23 

14.5108 

7 

14.2069 

24 

14.4164 

8 

18.0229 

25 

14.2250 

9 

17.9699 

26 

15.0203 

10 

17.9065 

27 

15.1082 

11 

17.7876 

28 

15.3061 

12 

17.5018 

29 

15.3798 

13 

16.8785 

30 

15.5854 

14 

14.4526 

31 

15.6785 

15 

14.2845 

32 

15.7742 

16 

14.2470 

33 

15.8130 

17 

14.2716 

34 

15.8130 

Data 

Pile:  D2209E 

Supply  Maaa  Plow  Rate 

(lba/aec) 

-  0.083142 

PORT 

PRBSSURB(pala) 

PORT 

PRESSURE (pa la) 

1 

114.0550 

18 

14.4190 

2 

114.9359 

19 

14.4074 

3 

19.3833 

20 

14.3492 

4 

19.3096 

21 

14.5936 

5 

19.2721 

22 

14.4345 

6 

14.2186 

23 

14.4669 

7 

14.2147 

24 

14.2949 

8 

19.2773 

25 

14.2263 

9 

19.2216 

26 

15.0707 

10 

19.1466 

27 

15.2182 

11 

19.0186 

28 

15.4910 

12 

18.6915 

29 

15.6164 

13 

17.9324 

30 

15.8906 

14 

14.4875 

31 

16.0302 

15 

14.2936 

32 

16.1970 

16 

14.2225 

33 

16.2436 

17 

14.2612 

34 

16.2578 

Data  Pile:  D2209G 


Supply  Mass  Plow  Rate 

(lbn/sec) 

«  0.108925 

PORT 

PRBSSURB(psia) 

PORT 

PRBSSURB(psia) 

1 

113.5972 

18 

14.4656 

2 

113.3916 

19 

14.4436 

3 

22.6264 

20 

14.3156 

4 

22.5579 

21 

14.6769 

5 

22.5605 

22 

14.3854 

6 

14.2225 

23 

14.3789 

7 

14.2108 

24 

13.9716 

8 

22.5346 

25 

14.2166 

9 

22.4777 

26 

15.1044 

10 

22.3950 

27 

15.3358 

11 

22.2462 

28 

15.8001 

12 

21.8583 

29 

16.1091 

13 

20.8510 

30 

16.6121 

14 

14.6169 

31 

16.9509 

15 

14.2897 

32 

17.3117 

16 

14.1487 

33 

17.4643 

17 

14.2095 

34 

17.5160 

Data 

Pile:  D2209H 

Supply  Mass  Flow  Rate  (lbm/sec)  *  0.12559 


PORT 

PRB8SURB ( ps ia ) 

PORT 

PRBSSURB(psia) 

1 

113.4084 

18 

14.4979 

2 

113.1343 

19 

14.4591 

3 

24.8079 

20 

14.3039 

4 

24.8661 

21 

14.7591 

5 

24.8726 

22 

14.3725 

6 

14.2056 

23 

14.3363 

7 

14.1992 

24 

13.8009 

8 

24.7187 

25 

14.2212 

9 

24.6424 

26 

15.0578 

10 

24.5558 

27 

15.3229 

11 

24.4161 

28 

15.9268 

12 

24.0359 

29 

16.4001 

13 

22.9484 

30 

17.0932 

14 

14.7268 

31 

17.5677 

15 

14.2781 

32 

18.0604 

16 

14.0893 

33 

18.2543 

17 

14.1888 

34 

18.2751 

Data  Pile:  D2209I 


Supply  Maes  Plow  Rate 

( lb*/ sec) 

-  0.13969 

PORT 

PRBSSURB(psla) 

PORT 

PRBS8URB<  psla 

1 

113.5312 

18 

14.5431 

2 

113.1923 

19 

14.4707 

3 

27.2235 

20 

14.2987 

4 

27.2597 

21 

14.8638 

5 

27.2700 

22 

14.3556 

6 

14.2082 

23 

14.3143 

7 

14.1992 

24 

13.6199 

8 

27.1446 

25 

14.2045 

9 

27.0683 

26 

14.9944 

10 

26.9726 

27 

15.2582 

11 

26.8135 

28 

15.9863 

12 

26.4256 

29 

16.5953 

13 

25.2256 

30 

17.5199 

14 

14.8716 

31 

18.1574 

15 

14.2677 

32 

18.8027 

16 

14.0169 

33 

19.0613 

17 

14.1733 

34 

19.1001 

Data 

Pile:  D2209J 

Supply  Mass  Plow  Rate 

(lba/sec) 

-  0.15997 

PORT 

PRBSSURB (psla) 

PORT 

PRBSSURB ( ps la 

1 

113.9075 

18 

14.5962 

2 

113.4640 

19 

14.4966 

3 

30.5765 

20 

14.2587 

4 

30.4433 

21 

14.9893 

5 

30.4562 

22 

14.3738 

6 

14.2056 

23 

14.2923 

7 

14.1953 

24 

13.4156 

8 

30.5274 

25 

14.2108 

9 

30.4343 

26 

14.9104 

10 

30.3140 

27 

15.0927 

11 

30.1343 

28 

15.9410 

12 

29.6778 

29 

16.7725 

13 

28.3161 

30 

18.0514 

14 

15.1613 

31 

18.9436 

15 

14.2755 

32 

19.8462 

16 

13.9315 

33 

20.1967 

17 

14.1552 

34 

20.2561 

Data  Pile:  DO 8 10 A 


Supply  Mass  Plow  Sate  ( lb*/ sec ) 


PORT 

PRBS8URB(psia) 

PORT 

1 

114.8473 

18 

2 

114.2306 

19 

3 

35.9310 

20 

4 

35.5198 

21 

5 

35.5172 

22 

6 

14.4174 

23 

7 

14.4161 

24 

8 

35.8961 

25 

9 

35.7939 

26 

10 

35.6491 

27 

11 

35.4564 

28 

12 

34.9469 

29 

13 

33.3719 

30 

14 

15.9291 

31 

15 

14.5118 

32 

16 

13.9752 

33 

17 

14.3527 

34 

Data 

Pile:  D0810B 

Supply  Maes  Plow  Rate 

(ltm/eec 

PORT 

PRBSSURB(psla) 

PORT 

1 

114.4594 

18 

2 

113.4785 

19 

3 

45.1444 

20 

4 

44.6841 

21 

5 

44.7410 

22 

6 

14.4122 

23 

7 

14.4135 

24 

8 

45.1302 

25 

9 

45.0035 

26 

10 

44.8289 

27 

11 

44.6026 

28 

12 

44.0492 

29 

13 

42.2272 

30 

14 

17.4950 

31 

15 

14.7769 

32 

16 

13.8859 

33 

17 

14.6902 

34 

0.19033 

PRBSSURB(psia) 

14.9023 

14.7652 

14.4859 

15.4467 

14.5002 

14.4717 

13.3545 

14.4161 

15.0717 

14.9838 

15.8605 

17.0062 

18.8877 

20.2571 

21.6433 

22.2304 

22.4102 


0.238907 

PRB8SURB  (pela) 

15.3601 

15.0859 

14.9398 

15.9252 

14.7226 

14.5726 

13.0014 

14.3941 

15.2269 

14.7368 

14.9489 

16.3299 

19.9067 

22.4270 

24.8878 

25.8886 

26.1330 


Appendix  B: 


This  appendix  Includes  samples  calculations  for  the  pipe 
calibration  and  for  the  first  step  of  the  numerical  simula¬ 
tion  code. 


Figure  B.l  Calibration  Manifold 

Pressure  Tap  1  -  Pressure  upstream  of  the  orifice  plate,  or 

pressure  upstream  of  flowmeter,  p* 

Pressure  Tap  2  -  Pressure  downstream  of  orifice  plate,  p* 

0  *  Mass  flow  rate 

g  *  Flowmeter  mass  flow  rate  reading 
A  «  surface  area  of  the  porous  pipe 
0  ■  porous  pipe  diameter  -  0.625  Inches 
Da.  -  Air  8upply  Pipe  Diameter  -  3.062  inches 
Da  ■  orifice  plate  bore  ■  1.0  inch 


L  ■  length  of  porous  pipe  saaple 
k  -  1.4 

C  ■  discharge  coefficient  ■  0.61 
3  -  Da/D».  -  0.326584 

R  «  ideal  gas  constant  for  air  ■  53.26  lbf-ft/lba-°R 

X  «  (px  “  Pa)/Px 

Y  -  1  -  (0.41  ♦  0.35  p*)x/k 

Saaple  Readings  froa  Orifice  Plate  Heasureaent  of  Hass  Plow 
Rate: 

Px  -  114.7224  psia 
Pa  ■  114.7007  psia 
pa  ■  25.19472  psia 
p«  -  25.16627  psia 
P*«»  *  14.3558  psia 
T  «  538  °R 
L  ■  1.953  inches 

X  -  0.000189 
Y  -  0.99994 

p  -  px/(RT)  -  0.5765  lba/ft’ 

Q  -  Da  C  Y  (  2  (px  -  Pa)/((1  -  «)J*'*/4 

Q  •  0.03604  lba/sec 
A  ■  no  L  *  0.02663  ft* 
pv  -  Q/(Ag.)  -  0.04206  lbf-sec/ft> 

AP*  ■  l(p»  ♦  p*)/2 J*  -  p.%»* 

AP*  ■  8,874,344.45  ibf*/ft« 


Maple  Readings  froa  Plowaeter  Mass  Plow  Rate  Heasureaent: 


px  -  15.00665  psia 
p»  ■  14.6B725  psia 
p«  -  14.68466  psia 
p«tt.  ■  14.37227  psia 
q  -  0.3975  cfa 
T  -  534  °R 
L  »  1.875  inches 

p  «  px/(RT)  «  0.07598  lba/ft* 

Q  ■  gp  -  0.0005034  lba/sec 
A  ■  no  L  -  0.025566  ft* 
pv  -  Q/(Ag.)  -  0.000612  lbf-sec/ft> 

APa  -  Up*  ♦  p« )/2 1*  -  p.*«*  -  189,011 . 3  lbf*/ft« 

Raaole  Calculations  for  Pirat  Increaent  of  Muaerical 

glaulatlon  Code  lx  -  Q.Q  to  x  »  ILiHSJUi 

(The  data  for  these  saaple  calculations  is  taken  froa  the 

raw  pressure  data  file  D2209B) 

p»  ■  pressure  at  the  upstreaa  end  of  the  evaporator  air 
supply  aanlfold 
pu  -  2791.20  lbf/ft* 

p«  •  pressure  at  the  downstreaa  end  of  the  evaporator  air 
supply  aanlfold 
p»  -  2777.88  lbf/ft* 

px  -  Internal  pressure  at  the  upstreaa  end  of  the  evaporator 
px  -  2781.648  lbf/ft* 
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pa  ■  pressure  at  the  first  pressure  tap  location  In  the 


evaporator 

p.  -  2775.931  lbf/ft* 

p>  ■  pressure  at  the  second  pressure  tap  location  in  the 
evaporator 

pa  «  2767.910  lbf/ft* 
k  -  1.4 

D  -  0.625  in  -  0.0521  ft 
PL  ■  total  pipe  length  -  6.976  ft 
A  *  constant  from  equation  2.1  ■  4.679  x  10*  ft*/sec* 

B  -  constant  fro*  equation  2.1  ■  4.2843  x  107  lbf/ft-sec 
T  -  533.1  °R 
x  «  0.005  ft 
Xi  -  0.00  ft 

xa  -  axial  location  of  first  pressure  tap  in  the  evaporator 
Xa  ■  0.50  ft 

Xa  ■  axial  location  of  second  pressure  tap  in  the  evaporator 
Xa  ■  0.99  ft 

Assuming  the  pressure  In  the  evaporator  changes  linearly  from 
upstream  value  to  its  downstream  value,  calculate  the  slope 
of  the  line. 

S  -  (p«  -  p«* ) / (PL/2 )  -  -3.8188  lbf/ft* 

Pi  »  Pi/(RT)  -  0.0978  lbm/ft* 

A.  «  n  D*/4  -  0.002132  ft* 

A»  »  no  ax  ■  0.0008164  ft* 

(Ap*)i  *  Pi*  -  P»l  ■  -53,251.844  lbf*/ft*  (negative  sign 


Indicates  mss  flow  into  pipe) 

Kx  -  C(X  -  Xa)(X  -  Xi ) 1 / I ( Xi  -  Xt)(Xi  -  Xa))  -  0.905 

Ka  -  Ux  -  Xi)(x  -  Xi)  )/( (Xa  -  Xa)  (X,  -  Xa)l  -  0.0201 

Ka  *  C (X  -  Xx)(X  -  Xa))/(<Xa  -  Xx)(Xa  -  Xa ) 1  «  -0.0051 

p(X  -  0.005)  -  Kip*  ♦  KaPa  ♦  K,p,  -  2781.6046  lbf/ft* 

Pm  -  p<x  -  0.005) 
p«.  ■  p.  4  8X  >  2791.18  lbf/ft* 

(Ap* ) a  -  P-2  -  pa.1  -  -53,365.507  lb£*/ft« 

(pv)  -  -q.l-B  ♦  (B*  ♦  4  A  |(Ap*)|  >*'*)/(2  A)  (fros  eqn  2.1) 
( pv ) x  •  -0.035672  lbB/ft*-sec 
( pv ) a  -  -0.035740  lbB/ft*-sec 
Max  ■  0.0 

u  -  axial  velocity 
ux  ■  0.0 

a*  -  q.kpx/px  «  1,281,139.47  ft*/sec* 
a  -  1,131.87  £t/sec 
Bx  ■  0.0  lbB/sec 

(Pv)-V,  -  ( (pv) x  ♦  ( pv) a  1/2  -  -0.035706  lba/ft*-sec 

pa  -  Pa/(RT)  -  0.0978  lbB/ft* 

Pmw  -  (Px  ♦  Pa )/2  -  0.0970  lbB/ft* 

Ua  ■  -(Pvl.vi  A»/(P.V«  Km)  ■  0.14015  £t/sec 

Ha,  -  ua/a  -  0.0001238 

u*«,  -  (ux  ♦  Ua )/2  «  0.07006  ft/sec 

Na.v,  -  (Max  +  Ma.)/2  -  0.0000619 

■  ■  *(pv).vt  A*  -  0.0000292  lbB/sec 
•  •  • 

m  -  Bx  ♦  ab  ■  0.0000292  IbB/sec 

•  e  • 

Bay,  ■  (Bx  ♦  Ba )/2  ■  0.0000146  ltMfc'eec 


Appendix  c:  Numerical  Simulation  code  *fid 
Sample  Results 


The  following  pages  are  a  listing  of  the  computer  code 
used  for  the  numerical  simulations.  The  code  is  written  in 
FORTRAN  77  and  some  sample  results  are  also  included  follow 
ing  the  program  listing. 


* 


Off B -DIMENSIONAL  HEAT  PIPS  SIMULATION  PROGRAM 


This  program  uses  one -dimensional  Incompressible  and  compres¬ 
sible  flew  models  to  solve  for  the  flow  properties  and  fric¬ 
tion  coefficients  In  a  closed  porous  pipe  with  blowing  and 
suction.  An  incompressible  model  consisting  of  Newton's  sec¬ 
ond  law  and  incompressible  flow  equations  is  used  to  calculate 
flow  properties  until  the  Mach  number  reaches  0.01.  After  that, 
the  compressible  model,  consisting  of  the  continuity  equation 
and  Shapiro's  method  of  influence  coefficients,  is  used  to 
solve  for  the  properties. 


DEFINITION  OF  VARIABLB8 : 

A  Constant  used  in  wall  boundary  condition  (ft*/sec*) 

AC  Cross-sectional  area  of  the  pipe  (ft*) 

AP  Surface  area  of  DX  increment  (ft*) 

B  Constant  used  in  wall  boundary  condition  (lbf/ft/sec) 

C2  Speed  of  sound  squared  (ft*/sec*) 

D  Pipe  inside  diameter  (ft) 

DMOM  Change  in  momentum  across  DX  increment  (lbf) 

DX  Special  step  size  (ft) 

DW  Change  in  mass  flow  rate  due  to  mass  transfer  at  the 
wall  (lbm/sec) 

F  Friction  factor 

F8  Shear  force  over  the  increment  DX  (lbf) 

FSSUM  Summation  of  all  the  incremental  shear  force  values 
over  the  entire  pipe  (lbf) 

FUD  Constant  to  account  for  the  statistical  variation  of 
the  porous  pipe  properties 
GAM  Ratio  of  specific  heats 

GC  Constant  in  Newton's  second  law  (lbm-ft/lbf/sec*) 

Ml  Mach  number  at  upstream  end  of  DX  increment 

M2  Mach  number  at  downstream  end  of  DX  increment 
M12  Ml  squared 

M2 2  M2  squared 

M8BAR  Average  of  M12  and  M22 

PATH  Atmospheric  pressure  (entered  in  psia  then  converted 
lbf /ft*) 

PCOND  Condenser  environment  pressure  (lbf /ft*) 

PBVAPD  Downstream  evaporator  environment  pressure  (entered  in 
psia  then  converted  to  lbf /ft*) 

PBVAPU  Upstream  evaporator  environment  pressure  (entered  in 
psia  then  converted  to  lbf /ft*) 

PBX  Environment  pressure  (lbf /ft*) 

PL  Pipe  length  (ft) 

P0  Upstream  end  of  porous  pipe  pressure  (lbf /ft*) 

PI  Static  pressure  at  upstream  end  of  DX  increment  (lbf /ft*) 


P2  Static  pressure  at  downstream  end  of  DX  Increment 
( lb£/f t* ) 

R  Ideal  gas  constant  £or  air  ( lbf-f t/lbm/*R ) 

RBW  Radial  Reynolds  number  “ 

REX  Axial  Reynolds  number 

RHO  Density  (lbm/ft*) 

RHOV1  Mass  flux  through  the  wall  at  the  upstream  end  of  the 
DX  increment  ( lbm/ftVsec ) 

RHOV2  Mass  flux  through  the  wall  at  the  downstream  end  of  the 
DX  Increment  ( lbm/f t*/sec ) 

RMU  Viscosity  ( lbf-sec/f t* ) 

8  Slope  of  the  evaporator  environment  pressure  line 

( lbf/f t*/f t ) 

TV  wall  shear  stress  (lbf/ft*) 

TO  Total  air  temperature  («R) 

VI  Mass  flow  rate  at  the  upstream  end  of  the  DX  lncresmnt 
(lbm/sec) 

W2  Mass  flow  rate  at  the  downstream  end  of  the  DX  incre¬ 
ment  (lbm/sec) 

XL  Local  axial  position  (ft) 

XND  Nondimens  Iona 1  axial  location  (XL/PL) 

X(l)-X(30)  Axial  location  of  pressure  measurements  (ft) 


DIMENSION  X( 30 ) ,P ( 30 ) 

REAL  K1,K2,K3,M1,H2,M12,M22 ,MSBAR 
X(l)  -  0.00 
X( 2 )  «=  6.0/12. 

X(3)  -  11.9/12. 

X(4)  -  18. /12. 

X( 5 )  -  23.9/12. 

X(6 )  «  30. /12. 

X( 7 )  -  35.9/12. 

X( 8  )  -  3.44633 
X(9 )  «  3.52967 
X(10)  «  43.84/12. 

X(ll)  »  45.78/12. 

X( 12 )  -  47.72/12. 

X( 13 )  -  49.72/12. 

X(14)  «  51.72/12. 

X( 15 )  -  53.66/12. 

X(16 )  «  55.66/12. 

X( 17 )  -  57.72/12. 

X(18)  -  59.72/12. 

X( 19 )  -  61.69/12. 

X( 20 )  -  63.72/12. 

X( 21 )  -  65.69/12. 

X( 22  )  -  67.66/12. 

X( 23 )  -  69.66/12. 

X(24)  -  71.66/12. 

X( 25)  -  73.66/12. 

X( 26 )  -  75.63/12. 

X( 27 )  -  77.66/12. 


X( 28 )  -  79.72/12. 

X(29 )  ■  81.72/12. 

X{  30 )  -  8.976 

OPEN ( 9 , PILE* • results ' ) 

OPBN( 10, FILE* 'data ' ) 

OPBK( 12 /PILE* ' fdat ' ) 

READ ( 10, 20 )  PATH 
20  FORMAT (P7 . 4  ) 

WRITE (9, 30)  PATH 

30  FORMAT (3X, 'PATH  -  ' ,F7 . 4, lx, 'psia • ) 
READ( 10,40)  TO 
40  FORMAT(F6 . 2 ) 

WRITB( 9 , 50 )  TO 

50  FORMAT ( 3X, 'TO  -  *,F6.2,'«R') 

RBAD( 10,60)  PBVAPU 
60  FORMAT(F7.4) 

WRITB(9, 70)  PEVAPU 

70  FORMAT(3X, 'PBVAPU  *  • ,F7 . 4, lx, • psia • ) 
READ ( 10, 80 )  PEVAPD 
80  FORMAT(F7 .4 ) 

WRITS (9, 90)  PBVAPD 

90  FORMAT (3X, 'PBVAPD  -  ' ,F7.4,lx, 'psia,/) 
DO  110  I  -  2,7 
RBAD( 10, 100 )  P(I) 

100  FORMAT (F7. 4) 

P(I)  *  P(I )*144. 

110  CONTINUE 

DO  115  I  -  10,30 
RBAD( 10, 112 )  P(I) 

112  FORMAT (F7. 4) 

P(I)  -  P( I ) *144 
115  CONTINUE 

WRITB(6, 120) 

120  FORMAT ( 3X, ' ENTER  THE  FUDGE  FACTOR') 
READ( 5, 130 )  FUD 
130  FORMAT(F5.3) 

WRITB(9, 140 )  FUD 


tv 

140  FORMAT(3X, 'FUD  -  ', 

m- 

PATH  -  PATM*144 . 

cm 

PBVAPU  -  PBVAPU*144 

PBVAPD  -  PEVAPD* 14 4 

PCOND  -  PATH 

V  -  ^ 

R  -  $3,335 

GAM  -1.4 

GAM1  •  (GAM  -  1 )/2 

D  «  0.625/12 

sew 

PL  -  6.976 

ft: 

PLH  -  PL/2 

>v 

DX  -  0.005 

ft. 

GC  -  32.174 

PI  -  3.1415927 

A  -  4679000000.0 

m: 

B  -  42843100.0 

FS6UM  «  0.0 
OW8UM  -  0.0 


ASSUMING  THE  EVAPORATOR  BMVI RON KENT  PRESSURE  CHANGES  LINEARLY 
PROM  ITS  UPSTREAM  VALUE  TO  ITS  DOWNSTREAM  VALUE,  CALCULATE  THE 
SLOPE  OP  THE  LINE. 


S  -  <  PBVAPD  -  PBVAPU ) /PLH 

CALCULATE  THE  PRBSSURB  IN  THE  PIPB  AT  THE  UPSTREAM  END  U8ING 
THE  PIR8T  THREE  MEASURED  PRBSSURB  VALUES. 


XI  -  X ( 2 ) 

X2  -  X( 3 ) 

X3  -  X ( 4 ) 

XC1  «  X2*X3/( (XI  -  X2 ) * (XI 
XC2  -  X1*X3/((X2  -  XI ) *( X2 
XC3  ■  X1*X2/((X3  -  XI ) * ( X3 
P(l)  -  P( 2 ) *XC1  ♦  P( 3)*XC2 
P0  *  P ( 1 )/144 


-  X3)) 

-  X3 )  ) 

-  X2)) 

♦  P ( 4 ) *XC3 


CALCULATE  THE  PRBSSURB  AROUND  THE  CONDENSER  ENTRANCE. 

P<8)  »  P( X*(PLH-0 . 5/12 ) ) ;  P{ 9 )  -  P(X-(PLH+0 . 5/12 ) ) 

XO  -  X( 9 ) 

XI  -  X(10) 

(•  X2  -  X<11) 

X3  «  X( 12 ) 

XC1  -  ((XO  -  X2)*(X0  -  X3) )/( (XI  -  X2)*(X1  -  X3)) 

XC2  »  ((XO  -  X1)*(X0  -  X 3 ) ) / ( (X2  -  X1)*(X2  -  X3)) 

XC3  -  ((XO  -  X1)*(X0  -  X2 ) )/( (X3  -  X1)*(X3  -  X2)) 

P(9)  -  P(10)*XC1  ♦  P(11)*XC2  ♦  P( 12 ) *XC3 
P(8)  «  P(9 ) 

DO  160  J  >  1,30 
PT  «  P(J)/144 
WRITE(9, ISO)  J,PT 

150  PORMAT( 3X, 'P( ' , 12, ' )  -  ' ,P7 . 4, IX, ' psl* ' ) 

160  CONTINUE 

WRITB(9,170) 

170  PORMAT(//,45X, 'RADIAL',9X, 'AXIAL* ,/,4X, 'X/L' ,7X, 'PRICTION' ,10X 
6 'MACH ' , 8X, 'REYNOLDS ' , 6X, 'REYNOLDS »,/, 15X, 'FACTOR ', 10X, 'NUMBER' 
CSX, 'NUMBER', 8X, 'NUMBER', 7X, 'BETA',/) 

WRITE( 12, 180 ) 

180  PORMATU2X,  'CALCULATED' ,9X,  '  BOWMANS  '  ,/,  4X,  *X/L',7X,  'PRICTION', 
C8X, 'FRICTION ',/,15x, 'FACTOR 10X, 'FACTOR' ) 

INITIALIZE  VARIABLES  AT  THE  UPSTREAM  PIPE  END 


XL  -  0.0 
W2  -  0.0 
M22  •  0.0 
P2  -  P(l) 
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T  -  TO 

P8X  «  PKVAPU 
RH02  «  P(1)/<R*T0) 

AC  -  PI*D**2/4 . 

AP  -  PI*D*DX 
I  -  1 
M  -  1 

C  -  AB8(P2**2  -  PBX**2) 

RHOV2  -  -PUD*GC* ( -B  +  SQRT(B*»2  ♦  4 .*A*C> )/(2.*A) 

C2  -  GC*GAM*P2/RH02 
XL  «  XL  ♦  DX 

BEGIN  MARCHING  DOVM  PIPE 

INCOMPRESSIBLE  MODEL  (M  <  0.01)  CALCULATIONS 

200  N  -  N  +  1 

IP( XL . LB . X ( 3 ) )  I  -  1 

I P ( XL . GT . X ( 3 ) . AMD . XL . LB . X ( 4 ) )  1-2 

IP(XL.GT.X( 1+2) )  I  -  I  +  1 

IP( XL . GB. X( 9 ) .AND. XL .LB . X( 11 ) )  1-9 

IP(XL.LT.X(8) )  PBX  -  S*XL  +  PBVAPU 

IP(XL.GT.X(9 ) )  PBX  -  PCOND 

W1  -  V2 

M12  «  M22 

Ml  -  SQRT(M12 ) 

U1  -  Ml*8QRT(C2 ) 

PI  -  P2 
RHOV1  -  RHOV2 
RHOl  ■  RH02 

IP( XL.GT.X( 8 ) . AMD. XL.LT.X( 9 ) )  GO  TO  210 

INTERPOLATE  TO  FIND  THE  PRB88URB  AT  8TATION  2 

XI  «  X(I) 

X2  -  X(I+1) 

X3  «  X( 1+2 ) 

CC1  -  P(I)/((X1  -  X2)*(X1  -  X3 ) ) 

CC2  -  P(I+1)/((X2  -  X1)*(X2  -  X3 ) ) 

CC3  -  P( 1+2 )/( (X3  -  X1)*(X3  -  X2 ) ) 

P2  -  CC1* ( XL-X2 ) * ( XL-X3 ) +CC2 * ( XL-X1 ) * ( XL-X3 ) +CC3* ( XL-Xl ) * ( XL-X2 ) 
XMD  «  XL/PL 
GO  TO  220 
210  P2  -  P(9 ) 

RHOV2  -  0.0 
XMD  -  XL/PL 
GO  TO  230 

CALCULATE  THE  PROPBRTXB8  AT  8 TAT I ON  2 

220  C  -  P2»*2  -  PEX**2 
ID  -  1 

IP(C.LT.O.O)  ID  -  -1 


C  -  ABS(C) 

RHOV2  «  ID*PUD*OC*(-B  ♦  SQRT(B**2  ♦  4*A*C) )/( 2*A) 

230  RHOVAV  -  (RMOV1  ♦  RHOV2 )/2 
RHG2  -  P2/(R*T) 

RHOBAR  -  (RHOl  ♦  RH02)/2 

U2  «  RH01*U1/RH02  -  RHOVAV*AP/(RH02*AC) 

M22  •  U2**2/C2 
M2  «  SORT (M2 2) 

UBAR  -  (VI  +  U2  )/2 
MSBAR  -  (M12  ♦  M22 )/2 
COM1  -  1  ♦  GAM1*M8BAR 
T  -  TO/COWl 
OV  -  -RHOVAV* AP 
DV8UM  =  DVB OH  +  DV 

IP( XL .GT. 3.489 . AMD . XL .LT. 3.491)  VCK  -  DW8UM 

V2  -  VI  ♦  DV 

VBAR  -  (VI  ♦  V2)/2 

RMU  «  2.27B-08*SQRT(T**3)/(T  ♦  198.6) 

RBX  -  4*VBAR/(PI*D*RMU*0C) 

RBV  «  RHOVAV*D/ ( RMU*GC ) 

BBTA  -  AB8 (RBV/RBX ) 

DMOM  -  V2*U2/GC  -  V1*U1/0C 

TV  «  ((PI  -  P2)*AC  -  DMOM )/AP 

PS  -  TV*AP 

P88UM  -  PS BUM  ♦  P8 

P  -  2*OC*TV/(RHOBAR*UBAR**2 ) 

ARBX  -  ABB (RBX) 

PSTAR  -  0 . 046/ARBX**0 . 2 
IP (XL. LB. 3. 488)  THBH 

PB  -  16/RBX*1 . 2337*BXP( 1 . 2* MSBAR) 

BL8B 

PB  -  PSTAR* (1+55*ARBX**.1*BXP(1.2*MSBAR)*BBTA**.9*PL/D**.l) 
BMD  IP 

CK  «  M/20  -  K 
IP(M.LT.IO)  00  TO  260 
IP(CX.BQ.O.O)  00  TO  260 
250  XL  *  XL  +  DX 

IP(XL.OT.PL)  GO  TO  390 

CHBCK  TO  8BB  IP  TUB  MACH  NUMBBR  HAS  RBACHBD  0.01.  IP  IT  HAS, 
SWITCH  TO  COMPRB88IBLB  PLOW  CALCULATIONS. 

IP (M2 .LB. 0.01)  OO  TO  200 
00  TO  300 

260  VRITS(9, 270)  XMD,P,M2, RBV, RBX, BBTA 

270  PORMAT(  3X,P5 . 3, 3X,P14 . 7, 5X,P9 . 7,-4X,P10 . 3, 4X,P9 . 2,4X,P7 . 4 ) 
VRITBt 12,282)  XMD,P,PB 
282  PORMAT( 3X,P5 . 3, 3X,P14 . 7,4X,r i r.B) 

IP(M.LT.IO)  00  TO  290 
K  •  K  ♦  1 


& 


OOMPRB8SIBLS  PLOW  CALCULATIONS 

M  -  M  4  1 

I P ( XL . LB . X ( 3 ) )  I  -  1 

IP( XL. OT.X( 3 ) . AMD. XL.LE. X(4 ) )  1-2 

IP(XL.OT.X( 1+2) )  I  -  I  4  1 

IP(XL.OB.X(9) . AND. XL. LB. X( 11) )  1-9 

IP(XL.LT.X(8) )  PBX  -  8*XL  4  PBVAPU 

IP(XL.GT.X(9) )  PBX  -  PCOMO 

VI  -  V2 

PI  -  P2 

RHOl  -  RH02 

RHOV1  -  RHOV2 

C2  -  0C*0AM*P1/RH01 

M12  -  H22 

Ml  «  SQRT(M12) 

U1  4  M1*8QRT(C2) 

IP ( XL .OT . X ( 8 ) . AMD.XL.LT. X( 9 ) )  00  TO  310 
INTERPOLATE  TO  PINO  THE  PRB88URB  AT  8TATI0N  2 


XI  -  X(I) 

X2  -  X(I41> 

X3  -  X(I42) 
PP1  -  P(I) 

PP2  -  P(I41) 
PP3  -  P(I42) 
XI  -  PP1/UX1 
K2  -  PP2/((X2 
K3  -  PP3/( (X3 


XI  -  PP1/UX1  -  X2)*(X1  -  X3)  ) 

X2  -  PP2/((X2  -  X1)*(X2  -  X3) ) 

X3  -  PP3/UX3  -  X1)*(X3  -  X2) ) 

P2  -  XI* ( XL-X2 ) * ( XL-X3 ) 4R2* ( XL-X1 ) • ( XL-X3 ) 4X3* ( XL-X1 ) * ( XL-X2 ) 
XMD  -  XL/PL 
00  TO  320 
310  P2  -  P(9) 

RHOV2  -0.0 
XMD  «  XL/PL 
00  TO  330 

CALCULATE  THE  PROPERTIES  AT  STATION  2 

320  C  -  P2**2  -  PBX** 2 
ID  -  1 

IP(C.LT.O.O)  ID  -  -1 
C  -  AB8(C) 

RH0V2  -  ID*PUD*OC* ( -B  4  8QRT(B**2  4  4*A*C) )/( 2*A) 

320  RHOVAV  -  (RH0V1  4  RHOV2)/2 
RH02  •  P2/(R*T) 

RHOBAR  -  (RHOl  4  RH02)/2 
DV  -  -RHOVAV* AP 
DVEUM  -  DVSUM  4  DV 

IP(XL.OT. 3. 489. AMD.XL.LT. 3. 491)  VCR  -  DVEUM 

V2  -  VI  4  DV 

VBAR  -  (VI  4  V2)/2 


V.'i'jV.'i* 


CM  -  M12*(1  ♦  GAM1*M12)*((W2/W1)**2)/((P2/P1)**2) 

M22  -  (-1  4  8QRT(1  4  4*GAM1*CM) )/< 2*0AM1 ) 

M2  «  8QRT(M22 ) 

MS BAR  -  (M12  4  M22 )/2 
C0W1  -  1  4  OAMl*MSBAR 
T  -  TO/ COW1 
C2  -  QC*GAM*P2/RH02 
U2  -  M2*8QRT(C2 ) 

UBAR  -  <01  4  U2>/2 

CFP  •  <1  4  GAM*M8BAR ) / ( GAM* MS BAR ) 

P  -  0 .5*  <  D/DX ) * ( ALOG(M12/H22 )  -  CPP*AL0G(P2/P1) ) 

TV  -  0 . 5*P*RHOBAR*UBAR**2/GC 

PS  -  TV*AP 

PSSUM  -  P88UM  4  P8 

RMU  «  2.27B-08*SQRT(T**3)/(T  4  198.6) 

RBX  «  4*WBAR/(PI*D*RMU*GC) 

ARBX  -  ABS(REX) 

RBV  ■  RHOVAV*D/ ( RMU*GC ) 

BBT A  -  AB8 ( RBW/RSX ) 

P8TAR  «  0.046/ ARBX* *0.2 
IP( XL .LB .3.488)  THBM 

PB  -  16/RBX*1 . 2337*BXP< 1 . 2* MS BAR ) 

BL8B 

PB  «  P8TAR*( 1+55* ARBX** . 1*BXP < 1 . 2*MSBAR) *BBTA** . 9*PL/D** . 1 ) 
BMD  IP 

CK  ■  M/20  —  K 
IP(CK.BQ.O.O)  GO  TO  360 
350  XL  -  XL  4  OX 

CHECK  TO  8BB  IP  THE  MACH  NUMBER  18  LBS8  THAW  0.01.  IP  IT  IS, 
SWITCH  BACK  TO  THB  I MCOMPRB8S 1 BLB  PLOW  CALCULATIONS. 

IP(M2.LB.0.01)  GO  TO  200 
IP (XL. LB. PL)  GO  TO  300 
GO  TO  390 

360  VRITS(9, 370 )  XMO,P,M2,RBV, RBX, BETA 

370  PORMAT( 3X,P5. 3, 3X,P14. 7, 5X,P9 . 7,4X,P10. 3, 4X,P9. 2, 4X,P7 .4) 
WRITB( 12, 382)  XMD,P,PB 
382  PORMAT( 3X,P5. 3, 3X,P14 . 7, 4X,P10. 8 ) 

K  ■  K  4  1 
GO  TO  350 

390  WRITB(9,400)  XWO,P,M2, RBV, RBX, BETA 

400  PORMAT( 3X,P5. 3, 3X,P14 . 7, 5X,P9 . 7, 4X,P10. 3, 4X,P9 . 2,4X,P7 . 4) 
VRITB( 12, 402 )  XWD,P,PB 
402  PORMAT( 3X, P5 . 3, 3X,P14 . 7, 4X,P10 . 8 ) 

COMPUTE  THB  PORCC  BALANCE  POR  THB  ENTIRE  TUBS. 

TEST  -  <P<1)  -  P<  30) )*AC  -  PSSUM 
WRITE (9,410 )  TE8T,OV8UM 

410  PORMAT(///, 3X, ' PORCE  BALANCE  - • ,P10.6, IX, * 1ME • ,/, 3X, 'MA88  BAL 


*,F9.6,2X, 'lba 


«ANCE  »' ,F6.4,lx, 'lba/sec* ) 

VRITE(9,  430)  VCX 

430  FORMAT ( 3X, *  Mass  Flow  Rate  Entering  Condenser  = 

a/sec ' ) 

CLOSE ( 9 ) 

CLOSE (10) 

CL08B( 12 ) 


TEST  RUM  II 


Supply  Mass  Flow  Rate  -  0.043651  lba/sec 


RADIAL 

AXIAL 

X/L 

FRICTIOM 

MACH 

REYNOLDS 

REYNOLDS 

BETA 

FACTOR 

NUMBER 

NUMBER 

NUMBER 

0.100 

0 . 099048 

0.018950 

-152.061 

7338.27 

0.0207 

0.200 

0.009700 

0.043203 

-200.811 

16726.01 

0.0120 

0.286 

0.014217 

0.071810 

-283.298 

27673.49 

0.0102 

0.372 

0.002657 

0.116745 

-454.632 

44442.84 

0.0102 

0.458 

0.002902 

0.191619 

-735.807 

71123.53 

0.0103 

0.515 

0.036181 

0.238237 

252.702 

86213.31 

0.0029 

0.630 

0.007313 

0.214018 

210.720 

76998.80 

0.0027 

0.716 

0.016780 

0.180273 

297.061 

65117.74 

0.0046 

0.816 

0.000037 

0.137835 

496.930 

50632.88 

0.0098 

0.917 

0.049359 

0.061729 

566.242 

18444.70 

0.0311 

Force 

Balance  ■  0 

.000287  lbf 

Mass  1 

Balance  ■  -0 

.0017  lba/sec 

Mass  1 

Plow  Rate  Inter ing  Condenser  ■  0.043609  lba^sec 

TIST 

RUM  12 

Supply  Mass  Flow 

Rate  «  0.05807  lbak/sec 

RADIAL 

AXIAL 

X/L 

FRICTION 

MACH 

REYNOLDS 

REYNOLDS 

BETA 

FACTOR 

NUMBER 

NUMBER 

NUMBER 

0.100 

0.141391 

0.018950 

-166.655 

7813.22 

0.0213 

0.200 

0.015890 

0.045054 

-238.355 

18550.60 

0.0128 

0.286 

0.012252 

0.078081 

-355.141 

31957.85 

0.0111 

0.372 

0.004333 

0.132968 

—602 .620 

53589.23 

0.0112 

0.458 

0.009795 

0.237445 

-1116.898 

91363.22 

0.0122 

0.515 

0.024286 

0.317003 

304.803 

115901.25 

0.0026 

0.630 

0.006685 

0.293726 

236.765 

106504.06 

0.0022 

0.716 

0.012464 

0.256543 

322.585 

93158.14 

0.0035 

0.816 

0.007668 

0.202693 

680.665 

75699.05 

0.0090 

0.917 

0.052697 

0.093339 

847.757 

35334.37 

0.0240 

Force  Balance  ■  0.000401  lbf 
Mass  'Balance  -  -0.0021  lbat/sec 

Mass  Flow  Rate  Inter ing  Condenser  *  0.058174  lba/sec 


* 


TB8T  RUM  93 


Supply  Mass  Flow  Rate  *  0.08314  lba/sec 


RADIAL 

AXIAL 

X/L 

FRICTIOH 

MACH 

REYNOLDS 

REYNOLD8 

BETA 

FACTOR 

NUMBER 

NUMBER 

NUMBER 

0.100 

0.257199 

0.017478 

-191.912 

8233.66 

0.0233 

0.200 

0.018961 

0.045229 

-298.903 

21284.82 

0.0140 

0.286 

0.011061 

0.082000 

-456.309 

38352.21 

0 .0119 

0.372 

0.007456 

0.147045 

-842.995 

67432.07 

0.0125 

0.458 

0.018558 

0.295485 

-1823.500 

124460.97 

0.0147 

0.515 

0.018951 

0.444090 

434.329 

167551.89 

0.0026 

0.630 

0.007233 

0.425127 

277.450 

158030.16 

0.0018 

0.716 

0.009127 

0.389829 

272.125 

143841.27 

0.0019 

0.816 

0.011444 

0.331650 

877.220 

127711.84 

0.0069 

0.917 

0.045913 

0.172545 

1303.856 

68633.47 

0.0190 

Force 

Balance  ■  0 

.001650  lbf 

Mass  1 

Balance  ■  0. 

0031  lb*/aec 

Mass  ! 

Plow  Rate  Entering  Condenser  -  0.082912  lbst/sec 

TE8T 

RUN  84 

Supply  Mass  Flow 

Rate  ■  0.108925  lbat/sec 

RADIAL 

AXIAL 

X/L 

FRICTION 

MACH 

REYNOLD8 

REYNOLDS 

BETA 

FACTOR 

NUMBER 

NUMBER 

NUMBER 

0.100 

0.277404 

0.015766 

-209.028 

0673.28 

0.0241 

0.200 

0.020153 

0.042822 

-353.672 

23576.27 

0.0150 

0.286 

0.009342 

0.081075 

-567.222 

44358.01 

0.0128 

0.372 

0.007461 

0.151141 

-1077.390 

81042.14 

0.0133 

0.458 

0.030118 

0.331755 

-2603.963 

157456.94 

0.0165 

0.515 

0.017825 

0.565882 

629.444 

223004.86 

0.0028 

0.630 

0.006598 

0.554932 

350.555 

213668.03 

0.0016 

0.716 

0.006495 

0.521664 

219.242 

197495.28 

0.0011 

0.816 

0.011422 

0.471743 

957.041 

187507.56 

0.0051 

0.917 

0.033950 

0.265377 

1754.761 

112193.95 

0.0156 

Force  Balance  -  0.024062  lbf 
Mass  Balance  -  0.0122  lha/sec 

Mass  Flow  Rate  Entering  Condenser  ■  0.108535  lb*/ sec 


TB8T  RUM  »5 


Supply  Hass  Plow  Rate  -  0.139692  1 be/sec 

RADIAL  AXIAL 


X/L 

PR  I CT I  OR 

MACH 

FACTOR 

MUMBBR 

0.100 

0.480220 

0.013170 

0.200 

0.000749 

0.044702 

0.372 

0.000610 

0.161377 

0.458 

0.041324 

0.362343 

0.515 

0.017733 

0.701334 

0.630 

0.006167 

0.706376 

0.716 

0.004396 

0.674508 

0.816 

0.009445 

0.635675 

0.917 

0.027973 

0.375817 

RBYMOLD8 

RBYMOLD8 

BBTA 

MUMBBR 

MUMBBR 

-274.625 

8659.52 

0.0317 

-500.176 

29494.15 

0.0170 

-1314.178 

104268.41 

0.0126 

-3468.318 

200480.73 

0.0173 

942.027 

294288.41 

0.0032 

439.304 

284733.22 

0.0015 

185.863 

264785.41 

0.0007 

896.208 

260769.89 

0.0034 

2200.172 

171322.94 

0.0128 

Pocce  Balance  ■  0.123370  lb£ 

Mass  Balance  «  0.0289  lb«/sec 

Hass  Plow  Rate  Inter Ing  Condenser  *  0.140039  ib«/aec 
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